Exact bounds on the ground-state energy of the infinite-f/ Hubbard model 
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We give upper and lower bounds for the ground-state energy of the infinite-f/ Hubbard model. 
In two dimensions, using these bounds we are able to rule out the possibility of phase separation 
between the undoped-insulating state and an hole-rich state. 
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The problem of phase separation (PS) in strongly cor- 
related systems is one of the most debated subject, espe- 
cially after the seminal paper by Emery and co-workers 
lyi, who also pointed out its relevance in the framework 
of high-temperature superconductivity. However, after 
several years of intensive numerical and analytical inves- 
tigation, there is no general consensus even for the quali- 
tative features of the possible instability in the Hubbard 
audi- J models ||-§. 

In this work we focus on the Hubbard model on a 
square lattice of L sites in the infinite-C/ limit, and prove 
that there is no PS in the low doping limit: our result be- 
ing one of the few rigorous statements in this field. The 
Hamiltonian is defined by 



W = -< E T^Gcl.Cj^aVG, 



(1) 



{hj),cr 



where the {i,j) denotes nearest neighbor sites, and the 
Gutzwiller projector Vg enforces the constraint of no 
double occupancy on each site. 

An upper bound on the ground-state (GS) energy 
can be easily obtained by use of a variational wave- 
function. A very simple choice is given by the fully 
polarized ferromagnetic (FM) state of energy EpM = 
2^ Sk X]i=i d *^*^^ ^«' where d is the dimensionality of the 
lattice and the sum over k is restricted to the lowest hole 
energy orbitals. 

A lower bound to the GS energy is instead more in- 
volved and can be derived as follows. The GS of the 
Hamiltonian Ti with M holes can be generally written as 



1^) 



J2 *(Ri...Rm)|Ri...Rm), 



(2) 



Ri...Ri\ 



where [Ri...Rjvf) are normalized states representing 
suitable superpositions of spin states with the con- 
straint of having the holes placed in the sites labeled by 
(Ri . . . Ra/). The sum is over all distinct hole configura- 
tions. The matrix elements Ktj of the hopping operator 
K on the reduced Hilbert space defined by the hole con- 
figurations are non vanishing only if a single hole hops at 



nearest neighbor distance. Therefore, Kij has no diag- 
onal elements. For instance, the element corresponding 
to the hopping Ri -^ R'l is (R'l . . . RM|i^|Ri . . • Rm) 
and its modulus is always smaller or equal than t. By 
construction, the lowest eigenvalue Eq of this matrix is 
the GS energy of the infinite-t/ Hubbard model. 

Now we prove that a lower bound of the GS energy 
is obtained by considering the matrix with all the non 
vanishing entries equal to — t. In fact, for the hermitian 
matrix with only non diagonal elements Kij , GS energy 
£^0 < and with components ipi = ^(Ri...Ra/) we 
have: 

>^B,j\^,\\t^j\>EB (3) 

Here, the matrix B describes a gas of M hard-core bosons 
(HCB) with nearest neighbor hopping Bij = —t for the 
non zero entries, which satisfies \Bij\ > \Kij\. The last 
inequality in (||) follows from the Perron-Frobenius theo- 
rem g and Eb is the lowest eigenvalue of B. This proves 
that a lower bound on the energy of the infinite-[/ Hub- 
bard model is given by the GS energy of an HCB gas 
with density equal to the doping of the infinite-[/ Hub- 
bard model: Eb < Eq < EpM- These bounds clearly 
hold in any dimension d. 

It is worth noting that the lower bound is not at all 
trivial, indeed i) the Hilbert space for a system of M 
HCB is much smaller than the Hilbert space of a system 
of spin-1/2 with M holes, ii) in general it is not true that, 
for a given Hamiltonian 7i, the energy of the bosonic 
GS is lower than the energy of the fermionic GS. This 
is actually the case for the spin-1/2 infinite-C/ Hubbard 
model but it does not hold in general (it is true if all 
the off-diagonal matrix elements of TL for the bosonic 
model are negative). For example let us consider the 
two-dimensional t — J model 
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where Si 


~ 2^(T,a' ^i,(j'''<y,<^'^i. 


matrices, 


and Ui = Yja 4,a'^h'y 



-fiin 



(4) 



, being t^.ct' the Pauli 

If we consider L — 16 
and 2 holes, the fermionic GS has an energy lower than 
the corresponding bosonic GS for J > 0.2t. Indeed for 
this Hamiltonian, even for bosons, the off-diagonal terms 
are not all negative definite. In Table I we report the GS 
energies of the t — J model for different value of J both 
for fermions and bosons. Remarkably Lanczos diagonal- 
izations on finite clusters show that, even at finite J, the 
hole-hole correlations in the t — J model are quite similar 
to the corresponding HCB results pOl]. 



-3.7 




In Fig. (P we show the previously discussed bounds 
on e{6). Since Eb{0) = EpMifi) = ^o(O) = 0, e(^) for 
the infinite-[/ Hubbard model lies in between tB{S) and 
^fm{&)- For any finite density the curvature of both 
curves is positive ruling out the possibility to have PS 
starting from 5 = Q. Notice that a previous lower bound 
to the GS energy given by Trugman |l3] just corresponds 
to e((5) = —4, and then it would not allow to draw any 
conclusion on the occurrence of PS in this model. 

In conclusion we have obtained a rigorous result that 
rules out PS in the infinite-L/ Hubbard model starting 
from the zero doping limit. Of course our findings do 
not exclude that PS could take place between two finite 
densities or at finite U > Q. 

This work has been partially supported by INFM and 
MURST (COFIN99). 



J/t 



£fh 



£b/t 



0.1 
0.2 
0.3 
0.4 
0.5 



-7.1348 
-8.2118 
-9.4199 
-10.6839 
-11.9843 



-7.5695 
-8.1967 
-9.3755 
-10.6056 
-11.8897 



TABLE I. GS energies £f and £b for 2 holes on L = 16 
sites for the fermionic and bosonic t — J model respectively. 



FIG. 1. e{5) for the FM state and for a gas of HCB. Full 
and empty circles are Monte Carlo data for L — 144 and 
L — 256, the dashed line is the analytical result for the FM 
state and the continuous lines are the low-density asymptotic 
behaviors. 

Let us consider now the relevance of these bounds on 
the GS energy in the two dimensional case of the infinite- 
U Hubbard model. The GS energy of an HCB gas can be 
computed numerically by quantum Monte Carlo on finite 
systems pi|. Moreover the T-matrix formalism in the 
low density limit allows to obtain the following analytic 
expression [p^: 



Eb{S) 
L 



-US - 27ri 



log 5' 



(5) 



with 5 = M/L. 

If PS occurs between a hole-free region (with S — 0) 
and a hole-rich phase (with S = 5c) then for finite size 
systems the function 



eiS) = 



Eo{S)-Eo{0) 
Mt 



(6) 



displays a minimum at J = (5c p| while, in the ther- 
modynamic limit, e(6) would be constant in the range 

0<6<Sr. 
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